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Pandemic H1N1/2009 viruses have been stabilized in swine herds, and some strains display higher 
pathogenicity than the human-origin isolates. In this study, high-throughput RNA sequencing (RNA-seq) is 
applied to explore the systemic transcriptome responses of the mouse lungs infected by swine (Jia6/10) and 
human (LN/09) H1N1/2009 viruses. The transcriptome data show that Jia6/10 activates stronger 
virus-sensing signals, such as the toll-like receptor, RIG-I like receptor and NOD-hke receptor signalings, as 
well as a stronger NF-kB and JAK-STAT singals, which play significant roles in inducing innate immunity. 
Most cytokines and interferon-stimulated genes show higher expression lever in Jia/06 infected groups. 
Meanwhile, virus Jia6/10 activates stronger production of reactive oxygen species, which might further 
promote higher mutation rate of the virus genome. Collectively, our data reveal that the swine-origin 
pandemic H1N1/2009 virus elicits a stronger innate immune reaction and pro-oxidation stimulation, which 
might relate closely to the increasing pathogenicity. 

A novel, reassorted H1N1 influenza virus outbreak began in early April 2009 (H1N1/2009) in Mexico and 
then rapidly spread worldwide in human populations, causing the first influenza pandemic of the 21 sl 
century 1 . Research reveals that hemagglutinin (HA ) and polymerase protein 2 (PB2) genes play important 
roles in the efficient transmission of the H1N1/2009 virus in animals 2 . The NA gene originated from Eurasian 
swine virus also contributes to the efficient respiratory droplet transmission of the H1N1/2009 virus through 
affecting the release of the viral particles 3 . Although the pandemic H1N1/2009 viruses display an efficient 
transmissibility that is equal to the seasonal influenza virus, the viruses are not characterized by a high mortality 
rate in infected people 4,5 . Subsequent laboratory researches also reveal that the patients infected by pandemic 
H1N1/2009 only develop non-complicated upper respiratory tract symptoms 5 . The data from animal model 
experiments, including studies in mice and ferrets, demonstrate that the pandemic H1N1/2009 viruses could 
stably replicate in the lung tissue but do not display high virulence in these animals 6 ". 

Analysis of inflammatory reaction in cell lines or in lung tissue infected by pandemic H1N1/2009 shows that 
although stronger inflammatory reactions are elicited by the H1N1/2009 virus than by the seasonal influenza 
virus, the reactions are much weaker than the highly pathogenic H5N1 avian influenza virus 9,10 . In the infected 
human macrophage, the pandemic H1N1/2009 virus does not elicit a "cytokine storm" 11 which is mainly 
responsible for pathological damage to the infected host cells. 

However, during the spread of pandemic viruses in human populations, viruses highly homologous with 
pandemic H1N1/2009 were isolated from swine herds in different countries, including Canada 12 , Thailand 13 , 
China 14 and other countries. These findings indicate that the pandemic H1N1/2009 virus could transmit into 
swine from people and has already stabilized in swine herds. 

In another alarming sign, the pandemic H1N1/2009 virus still undergoes further reassortment in swine. Virus 
containing HA gene originated from H1N1/2009 virus but neuraminidase (NA) and the internal genes originated 
from the European "avian-like" H1N1 and triple-reassortment swine H1N1 respectively was isolated from swine 
in Hong Kong in 2010 15 . Another reassorted virus containing the internal genes originated from the H1N1/2009 
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virus; however, the HA and NA genes originated from the swine 
influenza virus (SIV) (H1N2) was also isolated in pigs in the 
United Kingdom 16 . These isolates indicate that the pandemic 
H1N1/2009 viruses stabilized and underwent continuous reassort- 
ment in the swine herds. The novel viruses generated by the reassort- 
ment of H1N1/2009 with other influenza viruses in pigs will further 
alter their viral biological properties, affect the ability of swine-to- 
human or human-to-human transmission and potentially pose a 
significant threat to public health. 

In an epidemiological surveillance of SIV in central China carried 
out from August 2009 to April 2010, four viruses that are highly 
homologous with the H1N1/2009 virus were isolated from swine 
lungs. The four isolates caused mild symptoms with transient body 
temperature increases in the infected pigs. Meanwhile, two of them 
(Jia6/10 and 3/10) could cause severe weight loss and mortality in 
infected mice 14 . Virus A/Liaoning/ 1/2009 (H1N1) (LN/09) was iso- 
lated from the throat swabs of a 27-year old man who displayed mild 
flu-like symptoms. The sequenced genome of LN/09 (GenBank 
accession numbers: JX403975-82) displayed high homology with 
that of pandemic H1N1/2009 virus. In the present study, we showed 
that the swine-origin H1N1/2009 virus displayed higher pathogeni- 
city than the human-origin H1N1/2009 virus including more severe 
body weight loss and pathological lung damage as well as higher 
cytokine expression in infected mouse model. Given that the 
swine-origin pandemic H1N1/2009 viruses had displayed distinct 
biological properties and some of them displayed increased virulence 
as well as the potential for reverse swine-to-human transmission, we 
chose LN/09 and Jia6/10 to investigate the different pathogenesis and 



further to study the mechanism involving in the interaction between 
the human- and novel swine-origin H1N1/2009 isolates and their 
host. 

High-throughput RNA sequencing (RNA-seq) technology has the 
potential to reveal the alteration dynamics of the pathogen genome 
itself and the systemic change of host gene expression in the process 
of infection of pathogens, which could be beneficial to uncover the 
interaction mechanism and the pathogenesis of pathogens 1718 . In this 
study, by using Solexa/Illumina RNA-seq, transcriptome analysis 
was performed on mouse lung tissues that had been infected by 
Jia6/10 and LN/09 viruses. Our data revealed the distinct and subtle 
transcriptome landscape of mouse lung tissue infected by swine- and 
human-origin pandemic H1N1/2009 viruses. We also sought to 
obtain valuable information concerning the cellular gene expression 
patterns affecting lung function, which result in the different con- 
sequences for the two viruses. 

Results 

Increasing pathogenicity of the swine-origin H1N1/2009 virus. To 

investigate the differences of the pathogenicity of the swine- and 
human-origin H1N1/2009 viruses, the body weights and the 
pathological damage of the mice infected by Jia6/10 and LN/09 
were determined. As showing in Fig. 1A, the Jia/10 virus caused 
more significant body weight loss and four (50%) deaths of the 
infected mice. However, the body weight of the mice infected by 
LN/09 virus only displayed slight loss upto day 6 and then 
maintained increase. The virus titers of the Jia6/10 infected mouse 
lungs were higher than LN/09 infected mouse lungs on both day 3 
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Figure 1 | Pathogenicity of the Jia6/10 and LN/09 viruses in mouse model. (A) Body weight change of the Jia6/10- and LN/09- infected mice from day 0 
to day 10; (B) Virus titers of the Iia6/10- and LN/09-infected mouse lung tissues on day 3 and day 5; (C) Expression level of IL-ip, IL-6 and TNFain Jia6/ 
10 and LN/09 infected mouse lung tissues on day 3 and day 5. **p < 0.01, N represents the difference is not significant, and the p values are calculated 
using the student's t-test. 
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and day 5 (Fig. IB). Meanwhile, the mice infected by Jia6/10 deve- 
loped more severe pathological damage, including necrosis of the 
alveolar epithelial cells, edema and interstitium broadening as well 
as infiltration of the inflammatory cells. However, the pathological 
damages caused by the Jia6/10 and LN/10 viruses were much weaker 
than that caused by a highly pathogenic H5N1 influenza control 
(supplementary figure 1A [Supplementary Fig. S1A]). Further- 
more, the Jia6/10 virus stimulated higher expression of IL1|3, IL6 
and TNFa, indicating that this virus elicited a stronger innate immu- 
ne reaction in infected mouse lung tissues than LN/09 did (Fig. 1C). 

Analysis of the high-throughput sequencing (RNA-seq) tran- 
scriptome data. To gain a global and dynamic gene expression 
profile, RNA-seq was performed to explore the transcriptome of 
the mouse lung tissues infected by Jia6/10 and LN/09 on day 3 and 
day 5. More than 206 million 76-bp single-end reads were generated 
(Table 1). These valuable resources provide information to address 
two fundamental questions: Which level of the host transcriptome is 
affected by the virus infection, and what are the main differences in 
the transcriptome landscape among the samples infected by the two 
H1N1/2009 viruses? 

Generally, approximately 103 million reads could be mapped to 
the mouse genome, which represented approximately 50% of all the 
generated reads (Table 1). Of these reads, 94.8% mapped to anno- 
tated regions or predicted regions from the public annotation 
database, and 12.7% spanned exon-exon junctions. Notably, appro- 
ximately 13.6% of the reads were from pure intronic regions 
(Fig. 2A). Extensive reads mapping and assembling showed that 
although the RNA-seq data confirmed 67.2% and 65.6% of the anno- 
tated splice junction sites and exons, respectively, a considerable 
fraction of the detected splice junctions (15.6%) and the assembled 
exon blocks (24.5%) could not be incorporated into existing gene 
models (Fig. 2B). Most of these unincorporated junctions (77.7%) 
and exons (72.9%) overlapped with the annotated gene models, indi- 
cating the presence of many unidentified transcript isoforms or nat- 
ural antisense transcripts; whereas, the rest were mapped to 
intergenic regions, revealing many putative novel transcribed regions 
(NTRs) to be identified (Fig. 2B). 

To assess the relationship of the whole-transcriptome expres- 
sion patterns after the lung tissues were infected by the two viruses, 
a correlation matrix of pairwise comparisons was calculated 
(Fig. 2C). We found that the two points (day 3 and day 5) of the 
control and/or the Jia6/10-infected tissues were grouped separately. 
Notably, the LN/09-treated tissue on day 3 showed a higher correla- 
tion with the control group, while LN/09-treated tissue on day 5 was 
similar to the Jia6/10-infected group; however, the two points of the 
LN/09-treated tissues were also relatively closely matched (Fig. 2C). 
Meanwhile, some genome regions showed significant differences 
between the control and the virus-infected tissues, as exemplified 
by the local region of the Serum amyloid A-3 {SaaS) gene, which 
was greatly increased by virus infection. And the expression level 
of this gene was also different in Jia6/10 and LN/09 infected groups, 
especially on day 3 (Fig. 2D). 



Meanwhile, both viruses activated the production of numerous 
novel transcribed regions (NTRs) and alternative splicing (AS) 
events of genes. Notably, most of the NTRs were regulated by the 
virus infection and demonstrated diversity across the Jia6/10- and 
LN/09-infected groups (Supplementary Fig. S2; Table SI and 
Supplementary Materials and Results). In summary, a total of 337 
splicing events from 261 genes that changed significantly were iden- 
tified in the virus-infected groups (false discovery rate (FDR) < 0.05; 
Supplementary Table S2). Notably, a higher frequency of virus-spe- 
cific AS events was observed in the Jia6/10-infected tissue than in the 
LN/09-infected tissue on both day 3 and day 5 in all the eight differ- 
ent AS models (Supplementary Fig. S3A; Table S2 and Supple- 
mentary Materials and Results). Research has proved that a new 
splice variant of the human guanylate-binding protein 3 {GbpS) gene 
mediates anti- influenza activity through inhibiting viral transcrip- 
tion and replication 19 . And in this study, we identified the AS tran- 
scripts of Gbp5 gene in both Jia6/10 and LN/09 infected mouse lung 
tissue (Supplementary Fig. S3A). The detailed role of the AS of Gbp5 
gene in the pathogenesis of Jia6/10 and LN/09 infection need further 
investigation. Since AS plays an important function in inducing the 
biodiversity of transcriptome and proteome by producing multiple 
mRNA and protein isoforms with a single gene 20 , our results sug- 
gested that Jia6/10 exerted a greater influence on the transcriptome 
or proteome complexity than the LN/09 virus, and the complexity 
was increasing along the process of virus infection. 

Finally, further functional analysis of the genes involved in the 
virus-infected splicing events with the DAVID program (http://davi- 
d.abcc.ncifcrf.gov/) revealed that these genes were functionally 
enriched in different biological processes, including the alternative 
splicing, acetylation and phosphorylation. Interestingly, some of 
these genes participated in RNA splicing or RNA processing process, 
indicating that the virus promotes the RNA complexity in host cells 
by regulating the splicing regulator during their lifecycle. 

Dynamics of the virus-infected mouse transcriptome. To quantify 
and compare the gene expression level within and between the 
control and the virus-infected groups, we calculated and norma- 
lized the expression level of the genes through reads per kilobase of 
exon model per million mapped reads (RPKM). Using a threshold of 
> 1 RPKM in a given sample to define an "expressed" gene, 
2,246 genes (12.1%) showed highly differential expression between 
the control and the virus-infected mouse lungs on both day 3 and day 
5 with at least a fourfold change. Interestingly, the numbers of the 
upregulated genes were larger than the numbers of the downre- 
gulated genes in both virus-infected mouse groups at both time 
points (Fig. 3A; Supplementary Table S3). However, the distri- 
butions of the downregulated genes were much more diverse than 
the upregulated genes. Most upregulated genes were enriched in the 
common set of comparisons, while the downregulated genes were 
distributed within different single comparison groups (Fig. 3B). 
Consistent with this observation, the number of overall upregu- 
lated genes (1,060) was lower than that of the downregulated genes 
(1,153). Notably, approximately 100% of the differentially expressed 



Table 1 | Statistics of RNA 


seq data of the control and the Jia6/1 0 


and LN/09-infected groups 






sample 


Total reads number 


Mapped number (%) 


Junction num 


ber (%) * 


Control 3d (CK3d) 


34,768,246 


17,437,553 (50.15%) 


2,203,535 


12.64%) 


Control 5d (CK5d) 


36,749,042 


17,642,987 (48.01%) 


2,036,326 


1 1 .54%) 


Jia6/10 3d 


35,224,986 


14,497,494 (41.16%) 


1 ,844,470 


12.72%) 


Jia6/10 5d 


31,694,107 


19,422,629 (61.28%) 


2,461,232 


12.67%) 


LN/09 3d 


32,563,497 


16,735,295 (51.39%) 


2,193,341 


13.1 1%) 


LN/09 5d 


35,314,944 


19,601,355 (55.50%) 


2,338,884 


1 1 .93%) 


total 


206,314,822 


105,337,313 (51.06%) 


13,077,788 


; 12.42%) 


'Reads were aligned to the mouse genome by Bowtie and Tophat. 
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Figure 2 | Global overview of the RNA-seq data of the control and the Jia6/10- and LN/09-infected groups. (A) Pie chart demonstrating the 
proportion of RNA-seq reads assigned to mice annotated genomic features (left). The distribution of mapped reads from each sample is shown in the right 
bar chart. CK3d: control 3d group; CK5d: control 5d group; JIA6 3d: Jia6/10-infected 3d group; JIA6 5d: Jia6/10-infected 5d group; LN3d: LN/09-infected 
3d group; LN5d: LN/09-infected 3d group; the abbreviation of each group was also used in the subsequent figures. (B) RNA-seq data supporting for 
known and unknown gene models. Splicing junctions (left) are detected by TopHat and exons (right) are assembled by the Cufflinks program. The 
junctions and exons are then assigned to existing gene models in annotation or NTRs identified in this study. "Known" suggests that junctions or exons 
are identical with annotated gene features, while "unknown" means they are not incorporated into existing gene models. A diagram of gene model with 
two exons and one linked junction is shown in the middle. The numbers of annotated junctions and exons are indicated in parentheses. (C) Heat map 
matrix of pairwise Spearman correlation between six samples. (D) The distribution of RNA-seq reads (calculated as RPKM) and annotated genes 
(dashed red line) across the chromosome 7. An example of gene (Saa3) induced after virus-infection is shown on the bottom panel. 



genes (particularly for the upregulated genes) showed identical 
expression patterns across all investigated samples (Fig. 3C); with 
only 33 genes showed mixed expression patterns (Fig. 3D, inserted 
pie chart). Furthermore, to examine the biological roles of the 
differentially expressed genes, a Gene Ontology (GO) enrichment 
analysis was applied to the up- and downregulated genes. Notably, 
the upregulated genes were mostly enriched in the regulation of 
immune reactions, such as the immunity system progress, the im- 
mune response, the defense response, the innate immune response 
and the inflammatory response. However, most downregulated 
genes were involved in the cellular metabolism processes, such as 
striated muscle contraction, oxygen transport and the regulation of 
biological reaction (Fig. 3D; Supplementary Table S4). The 



downregulated genes engaged in the regulation of oxygen 
transport and gas transport indicated that the respiratory and gas 
metabolism function weakened or were damaged after virus 
infection, which might be closed related to the pathogenesis of the 
virus. To directly determine which pathways the two viruses affected 
after infection, all of the differentially expressed genes were mapped 
to the entire pathway based on the Reactome database 21 . The 
upregulated genes were observed largely enriched in the highly 
connected or pivot nodes in the pathway, while the downregulated 
genes were distributed at the terminal nodes of gene network. Two of 
the largest clusters of upregulated gene were found to participate in 
cell cycle regulation and interferon (IFN) regulation pathways 
(Fig.3E). 
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Figure 3 | Analysis of the differential genes and the functional maps of the Jia6/10 and LN/09-infected transcriptome. (A) Differentially expressed 
genes (with at least a 4-fold change) between control (CK) and virus-infected cells. The number of up- (pink) or down-regulated (green) genes for each 
comparison are shown in the bar chart. (B) Venn diagram showing the distribution of shared differentially expressed genes (indicated in (A)) among 
JIA63d, JIA65day, LN3d and LN5d. Gene number is listed in each of the categories. Numbers in parentheses indicate the number of up- or 
down-regulated gene within each category. Please note that the total number of up- and down-regulated genes in the parenthesis may not be equal to the 
number of genes in some categories, because some genes showed mixed expression patterns (i.e., up-regulated in one comparison, while down-regulated 
in another comparison, and vice versa) in the categories. (C) Two-way clustering analysis of all the 2,246 differentially expressed genes. (D) Gene 
Ontology (GO) enrichment analysis for the up- (top bar chart) and down-regulated genes (bottom bar chart), respectively. Only the top 15 enrichment of 
GO terms from "biological process" category are listed here. The inserted pie diagram depicts the number of up- and down-regulated or mixed genes. 
(E) Pathway analysis of differentially expressed genes. All the differentially expressed genes are mapped to the entire pathway based on Reactome. Inserted 
boxes show the two largest sub-pathways enriched within the up-regulated genes. 



Two virus infections impact the cellular innate immunity re- 
sponse through distinct models and with different magnitudes. 

As mentioned above, although the change in magnitude was not 
identical, the upregulated genes triggered by the Jia6/10 and the 
LN/09 infections displayed identical change patterns and were asso- 
ciated with biological functions related to innate immune processes. 
Therefore, the differences between the innate immune reaction 
caused by the Jia6/10 and LN/09 infection were comprehensively 
discussed. 

The innate immune system is the first defense of an organism and 
is associated with a highly conserved host-cell signaling mechanism 
that protects the host against invading pathogens 22 . Different pat- 
tern-recognition receptors (PRRs) are expressed to recognize patho- 
gen-associated molecular patterns (PAMPs), which initiate the 
signaling cascades inducing the production of proinflammatory 
cytokines and IFNs 22 . Currently, three types of PRRs have been iden- 
tified (to be activated in response to viral pathogens infection): the 
RIG-I like receptors (RLRs), the toll-like receptors (TLRs) and the 



NOD-like receptors (NLRs). The expression levels of all TLR genes 
(except TlrS) increased in response to Jia6/10 and LN/09 infection 
(Fig. 4; Supplementary Table S3). Among the different TLRs, TLR3 
and TLR7 mainly recognize dsRNA and ssRNA, respectively, which 
would be produced during influenza virus infection 23,24 . In this study, 
both Jia6/ 1 0 and LN/09 activated the TLR signaling mainly through a 
Tlr7/Myd88/Irf7-dependent mechanism. Previous study indicated 
that mice lacking Tlr7 and Myd88 fail to induce antiviral protection 
against influenza virus challenge 25 . Our data showed that the degree 
of increase of these genes was apparently higher in the lungs of the 
Jia6/10-infected mice than in the LN/09-infected groups on both day 
3 and day 5. The magnitude of the signal activation decreased from 
day 3 to day 5 in both the Jia6/10- and LN/09-infected mouse groups 
(Fig. 4). The RLR family, which includes RIG-I, Mda5 and Lgp2. 
RIG-I and Mda5 could recognize single-stranded 5 '-triphosphate 
RNA with base-paired structures and long dsRNA, respectively 26,27 . 
The RLR and the TLR signalings sensing influenza virus infection 
play partially redundant roles in inducing innate immunity 28 . In the 
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Figure 4 | Comprehensive analysis of the gene expression pattern of the innate immune and oxidation reaction after the Jia6/10 and the LN/09 
infection. Genes play significant role in the TLR, RLR and NLR signalings as well as NF-kB and JAK-STAT signalings are included (upper and middle 
left) . The up-regulated genes which encode interleukins, interferons, chemokines and ISGs as well as those associated with influenza virus infection are 
analyzed (lower left). Genes associated with pro-and anti-oxidation reaction also listed in the figure (right part). The six grids in each line represent the six 
groups (CK3d, CK5d, JIA63d, JIA6/5d, LN3d and LN5d), and the color in the grid represents the relative expression levels showing in the lower left 
quarter. 



present study, both RIG-I and Mda5 genes were upregulated in the 
Jia6/10- and LN/09-infected mouse lungs on day 3 and day 5. The 
expression level of the two genes was higher in the Jia6/10 groups 
than in the LN/09 groups. Furthermore, the expression level of the 
two molecules decreased in the Jia6/10-infected mouse lungs but 
maintained an increased level in the LN/09-infected groups from 
day 3 to day 5, indicating that the activation of RIG-I and Mda5 
did not progress at the same rate during the Jia6/10 and the LN/09 
infections (Fig. 4; Supplementary Table S3). Meanwhile, a signifi- 
cant activation of Lgp2, a negative regulatory gene of RIG-I and 
Mda5 29 , was observed in both the Jia6/10- and the LN/09-infected 
mouse lungs (Fig. 4). The continuous increase of Lgp2 from day 3 to 
day 5 in both the Jia6/10- and the LN/09-infected groups indicated 
that the organism prevented the excessive production of the down- 
stream inflammatory cytokines by regulating the activation of the 
RLR signaling. NLRP3, an NLR signal, could form a functional 
inflammasome complex with PYCARD (its adaptor protein), which 



responds to external stimuli and plays a role in regulating the activa- 
tion of capase-1 and the maturation of IL-lb and IL18 30 . In the 
present study, four members of the NLR signaling group {Nod2, 
Nlrp3, Pyard and caspase 1 genes) increased after the Jia6/10 and 
the LN/09 infections (Fig. 4). Generally, the Jia6/10 virus activated 
stronger TLR, RLR and NLR signaling reactions than LN/09, which 
might be associated with the higher pathogenicity of Jia6/10. 

NF-kB is a protein complex that controls the transcription of a 
large number of genes associated with inflammatory cytokines and 
the downstream interferon-stimulated genes (ISGs) 31 . The express- 
ion of Nfkbl, Ikbka, Ikbkb and Ikbkg genes in the canonical pathway 
of NF-kB signaling showed minor changes after infection with both 
viruses. This result was consistent with the upregulation of the NF- 
kB inhibitors {Nfkbia and Nfkbib) after the Jia6/10 and the LN/09 
infections (Fig. 4). However, the expression of Nfkbl and Ikbke genes 
in the non-canonical pathway of NF-kB signaling increased in both 
virus-infected groups (Fig. 4), which might suggest that the Jia6/10 
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and the LN/09 viruses both activated NF-kB signaling, with different 
magnitudes, through the non-canonical pathway but not through the 
canonical pathway. 

The cytokines and chemokines play significant roles in defending 
against the virus infection and recruiting inflammatory cells to the 
infection site. The data revealed that the important pro-inflammat- 
ory factor genes (ILlb and IL6) increased significantly in both virus- 
infected groups, and the expression values were obviously higher on 
day 3 than day 5. However, the expressions levels of most chemo- 
kines maintained an increase from day 3 to day 5 (Fig. 4; Supple- 
mentary Table S3). The different trends might relate to the different 
production processes associated with the pro-inflammatory cyto- 
kines and the chemokines after virus infection. Notably, certain che- 
mokines (e.g., CxcllO) significantly increased in response to the Jia6/ 
10 (294-fold on day 3 and 537-fold on day 5) and LN/09 (57-fold on 
day 3 and 639-fold on day 5) infections. Meanwhile, chemokines 
Ccl2, Ccl3, Ccl4, Ccl5, Ccl7, Cxcll, Cxcl2, Cxcl9, CxcllO, Cxcll3 and 
Cxcll6 were more highly expressed in the Jia6/10-infected groups. 
The IFNs, e.g., Ifnb, Ifng, 1128a and 1128b, as well as tumor necrosis 
factor alpha {Tnfa), were significantly upregulated in the two virus - 
infected mouse lungs, particularly the Ifng in the Jia6/10 (215-fold) 
and the LN/09 (152-fold) infection groups on day 5. Notably, the 
expression of most cytokines and chemokines was higher in the Jia6/ 
10-infected mouse lungs than in the LN/09-infected groups, which 
indicated that Jia6/10 activated a much stronger inflammatory res- 
ponse and greater inflammatory cell recruitment (Fig. 4; Supple- 
mentary Table S3). 

By binding to the receptor, type I IFN will activate the JAK-STAT 
signaling pathway and promote the expression of numerous ISGs, 
which play a significant role in defending against virus infection 32 . In 
this study, the data showed that the expression of Jak2, Statl and 
Stat2 genes increased in response to the Jia6/10 and the LN/09 infec- 
tions. The expression of different types of ISGs, such as Mxl, Isgl5 
and members of the Oas gene family also increased after infection by 
Jia6/10 and LN/09. Notably, the expression level of most ISGs was 
higher on day 3 than on day 5 in both the Jia6/10- and LN/09- 
infected mouse lungs. The expression level was obviously higher in 
the Jia6/10-infected groups than in the LN/09 groups (Fig. 4; 
Supplementary Table S3). This difference might be associated with 
the virus replication process and the higher virulence of Jia6/10 in 
infected mice. 

Pro- and anti-oxidation effects elicited by virus infection. The 

production of reactive oxygen species (ROSs) and reactive nitrogen 
species (RNSs), particularly superoxide, is an important host defense 
mechanism for clearing invading pathogens 33 . However, excessive 
oxidation stress is an important consequence of the immunopath- 
ology associated with influenza virus infection 34 . The important pro- 
oxidation genes, Cybb, Xdh, Nos2 and Ncfl, which are responsible for 
generating free radicals, increased in the Jia6/10- and the LN/09- 
infected mouse lungs on day 3 and day 5. Meanwhile, most anti- 
oxidation genes, such as Sodl, Nqol, Cat, Dhdh, Gstpl, Gstp2 and 
Gpx4, decreased in both the Jia6/10- and the LN/09-infected mouse 
lungs (Fig. 4; Supplementary Table S3). Overall, the increase in the 
pro-oxidation factors and the decrease in the anti-oxidation factors 
were more significant in the Jia6/10-infected groups than in the LN/ 
09 groups, indicating that the Jia6/10 virus stimulated greater ROS 
and RNS production, which might be associated with its higher 
virulence in mice. 

Meanwhile, the ROS-generated DNA damage presents a signifi- 
cant mutation risk and might accelerate organism evolution 35,36 . In 
this study, we also mapped the sequencing reads to the genomes of 
Jia6/10 and LN/09. The results revealed that a total of 29 point 
mutations occurred in the Jia6/10 virus genome, while the mutation 
number in LN/09 virus genome was only 19. Six mutations occurred 
in the HA gene of the Jia6/10 virus; however, no mutation was found 



in the LN/09 HA gene (Supplementary Table S5). Notably, the 6 
mutated nucleotides resulted in 5 amino acid mutations, and 2 
located at the antigenic site of the HA protein. The high mutation 
rate of the Jia6/10 genome might relate to the stronger pro-oxidation 
reaction elicited by the Jia6/10 virus and may partially verify the role 
of ROS in the promotion of genome mutation. 

Discussion 

Previous studies have shown that the pandemic H1N1/2009 stabi- 
lized in the swine population in different countries 1214 . The human- 
to-swine transmission of the pandemic H1N1/2009 virus could 
substantially affect the evolution and epidemiology as well as the 
biological properties of the virus. In this study, we have showed that 
some swine-originated H1N1/2009 viruses displayed increasing 
pathogenicity compared with the human H1N1/2009 virus. Jia6/10 
could cause mortality and significant weight loss of the infected mice 
without previous adaptation, while the LN/09 virus only raises the 
transient weight loss after infection. Meanwhile, the Jia6/10 virus 
results in more severe pathological damage (compared to LN/09) 
in infected mouse lung tissues. However, the mechanism for the 
increased pathogenicity after the transmission of pandemic H1N1/ 
2009 to swine and the interaction between swine pandemic H1N1/ 
2009 and the host is still unclear and deserves clarification. In an 
attempt to reveal the mechanism underlying the increasing patho- 
genicity after the transmission of pandemic H1N1/2009 to swine, 
RNA-seq was performed to explore the systemic and subtle tran- 
scriptome alterations of the mouse lung tissues infected by the swine 
and human pandemic H1N1/2009 viruses. 

The global transcriptome data showed that the gene expression 
patterns, for both protein-coding genes and newly identified NTRs, 
significantly changed after the lung tissue was infected by Jia6/ 1 0 and 
LN/09 viruses. Meanwhile, numerous NTRs and ASs might contrib- 
ute to the extensive level of transcriptome complexity due to virus 
infection. Notably, the differential upregulated genes were mainly 
associated with immune reaction regulation, and the activated innate 
immunity might be associated with pathological injury to the lung 
tissue. 

Acute lung injury (ALI) is a common consequence of a cytokine 
storm in the lung alveolar environment and systemic circulation 37 . 
The balance of pro- and anti-inflammatory mechanisms is critical to 
maintaining lung immune homeostasis. In the present study, the 
three innate immunity regulating pathways (TLR, RLR and NLR 
signalings) and numerous inflammatory cytokines were activated 
in response to the infection of Jia6/10 and LN/09, particularly on 
day 3 in the Jia6/10-infected group. All these data showed that the 
Jia6/10 virus activated a much stronger inflammatory response than 
LN/09 on both day 3 and day 5. The changing trend of the activated 
virus-sensing signaling is consistent with the virus titers of Jia6/10- 
and LN/09 infected mouse lung tissues and the pathogenicity of the 
two viruses: the higher virulent Jia/10 virus stimulated stronger 
innate immune signaling on day 3 which in turn inhibited its further 
replication on day 5, while the lower virulent LN/09 virus stimulated 
a relatively weak innate immune signaling on day 3 which nearly did 
not affect its further replication on day 5. 

However, when compared with the transcriptome data from the 
lungs of mice infected with a highly pathogenic H5N1 influenza 
virus, (A/duck/Hubei/hangmeiO 1/2006 [Hm/06]), which could lead 
to a cytokine storm and result in severe pathological damage to the 
lungs 38,39 , the expression levels of most chemokines in the Jia6/10 and 
LN/09 groups were obviously lower than in Hm/06 groups (data not 
show). Meanwhile, the pathological damage associated with Jia6/10 
infection was obviously weaker than that of Hm/06 (Fig. 1), which 
might indicate that both Jia6/10 and LN/09 didn't promote a cyto- 
kine storm in the lungs of infected mice. This observation is consist- 
ent with the upregulation of distinctive negative regulatory genes, 
such as Tollip, Socsl and Socs3, which involve in inhibiting different 
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cytokine-induced signaling. Meanwhile, the activated JAK-STAT 
signaling, as well as the upregulation of ISGs, might play a significant 
role in restricting virus replication. This indicates that the pathogen- 
infected tissues are in a state of "high-alert" situation in which the 
inflammatory reaction can either be efficiently initiated if the infec- 
tion progresses or prevented if the infection resolves 40 . 

Apart from stimulating a robust pro-inflammatory response, the 
virus infection also results in the production of ROSs (0 2 ~, OH~, 
H 2 0 2 ) and RNS intermediates (nitric oxide (NO), N0 2 , HN0 2 , 
ONOO~) 41 . While both ROSs and nitrogen intermediates display 
potent antimicrobial activity against a variety of pathogens, they 
can also cause extensive tissue damage if produced in excess 34 ' 41 . 
Upregulated pro-oxidation genes and downregulated anti- oxidation 
genes might elicit excessive amounts of ROS and RNS in both Jia6/10 
and LN/09 infections. Although the pro-oxidation ability is relatively 
weaker than the highly pathogenic H5N1 virus hm/06 (data not 
show), the Jia6/10 virus elicits a much stronger pro-oxidation 
reaction than LN/09. Previous study has verified that the oxidized 
phospholipid (OxPL) was produced through TLR4-TRIF-TRAF6 
pathway and was able to induce ALI in inactivated H5N1 challenged 
mice. Tlr4 silencing could protect mice from H5Nl-induced ALI, 
and inhibiting the expression of Ncfl , which regulates the production 
of ROS, relieves the severity of H5Nl-mediated ALI 42 . The results are 
well consistent with our present observation that the mice infected by 
Jia6/10 displayed more severe lung damage accompany with higher 
expression of Tlr4 and Ncfl in this study. 

The excessive production of ROS and RNS is harmful to organ- 
isms, whereas, they might be beneficial for the mutation of virus 
genomes, which could produce more virulent and adaptable viruses. 
This phenomenon might be one of the potential causes of the 
increasing virulence of the Jia6/10 virus. 

A virulent H1N1/2009 virus (MA-CA/04), containing 5 amino 
acids mutations (G155E, S183P and D222G in HA; E158G in PB2; 
D101G in NP) was selected through serial lung passage of A/CA/04/ 
2009 (CA/04) in mice 43 . However, none of these sites was found in 
the mutations of Jia6/10 and LN/09 genome. Interestingly, the ori- 
ginal amino acids at site 222 of the HA protein of Jia6/10 and LN/09 
viruses were G and D, respectively, indicating the potential high 
virulence of Jia6/10 virus. Notably, we identified the G222E mutation 
in the genome of Jia6/10 infected mouse lung tissues, suggesting that 
the mutation of the virus genome is complex and multi-directional. 
Further microarray analysis of the mouse lung tissue infected by CA/ 
04 and MA-CA/04 revealed that MA-CA/04 infection was associated 
with stronger immune response and some of the differential genes 
enriched in interferon regulatory factor (IRF) binding motif, 
immune cell trafficking genes including chemokines and their recep- 
tor as well as proinflammatory cytokine genes such as IL6 and Tnfa 4i . 
Most of the differential genes associated with immune response were 
also identified in Jia6/10 infected mouse lung tissues in the present 
study. The similar expression pattern of these genes related to 
immune reaction suggests their important functions in the increas- 
ing pathogenicity of the Jia6/10 and MA-CA/04 viruses. 

In conclusion, both Jia6/10 and LN/09 elicit a complex cellular 
response in infected mouse lung tissues. The innate immune res- 
ponse and the pro-oxidation process are activated more significantly 
on day 3 rather than on day 5 in mouse lungs infected with the swine- 
origin pandemic H1N1/2009 virus. Although the inflammatory reac- 
tion and pro-oxidation are still weaker than those produced by the 
highly pathogenic H5N1 avian influenza virus and neither virus 
elicits a cytokine storm, Jia6/10 triggers a stronger inflammatory 
reaction and pro-oxidation stimulation and higher cytokine produc- 
tion than LN/09. These differences might be closely related to the 
increasing virulence of Jia6/10. In light of the increase in virulence 
and the potential for reassortment and reverse swine-to-human 
transmission, it is essential to investigate the epidemiology and the 
evolution of the swine-origin pandemic H1N1/2009 virus. 



Methods 

Pathogenicity of the swine- and human-origin H1N1/2009 viruses. Stock viruses 
were prepared by propagating the LN/09 and Jia6/10 viruses in MDCK cells at 37°C in 
a 5% C0 2 incubator. The infectious titers were determined using the Reed and 
Muench method 44 . Each mouse was infected intranasally with 50 uL of virus stock 
fluid diluted in PBS to 10 3 TCID 50 . The change of the body weight and the 
histopathological analysis of the infected mouse lung tissues on day 3 and day 5 were 
determined after Jia6/10 and LN/09 infection. The production of cytokines IL-lp\ 
IL-6 and TNFa were also detected by enzyme-linked immunosorbent assay (ELISA). 
The detailed procedures were described in Supplementary Materials and Results. 

All animal experiments were conducted in compliance with Guide for the 
Management and Use of Laboratory Animals (Hubei province, July, 2005) and the 
legislation of Experimental Animal Management Ordinance (Ministry of S&T of the 
people's Republic of China, Oct, 1988). The animal studies were approved by 
Experimental Animal Center, Institute of Medicine, Hubei Province, China (No. 
00018734). 

RNA preparation and Solexa/Illumina sequencing. Each mouse was infected 
intranasally with 50 uL of Jia6/ 1 0 and LN/09 viruses at 1 0 3 TCID 50 . On day 3 and day 
5, the infected- and control- mouse lung tissues were collected and applied to RNA 
extraction and Solexa/Illumina sequencing. The detailed procedures and the Solexa/ 
Illumina sequencing were described in Supplementary Materials and Results on line. 

Bioinformatics analysis. To systematically characterize the transcriptome of virus- 
infected tissue, we employed various publicly available tools for the mapping, 
assembly and quantification of transcripts. We also integrated these tools with 
additional informatics filtering steps to enrich the results for the most robust 
transcriptome construction. Read alignment and assembly, discovery of splice 
junctions, identification of NTRs and AS events as well as the quantification of gene 
expression levels were described in Supplementary Materials and Results on line. 
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